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Abstract 

The cooling of neutron stars by URCA processes in the kaon-condensed neutron 
star matter for various forms of nuclear symmetry energy is investigated. The kaon- 
nucleon interactions are described by a chiral lagrangian. Nuclear matter energy is 
parametrized in terms of the isoscalar contribution and the nuclear symmetry energy 
in the isovector sector. High density behaviour of nuclear symmetry energy plays 
an essential role in determining the composition of the kaon-condensed neutron star 
matter which in turn affects the cooling properties. We find that the symmetry 
energy which decreases at higher densities makes the kaon-condensed neutron star 
matter fully protonized. This effect inhibits strongly direct URCA processes result- 
ing in slower cooling of neutron stars as only kaon- induced URCA cycles are present. 
In contrast, for increasing symmetry energy direct URCA processes are allowed in 
the almost whole density range where the kaon condensation exists. 
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star cooling 
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1 Introduction 



The possibility of the kaon condensation in dense nuclear matter, proposed 
by Kaplan and Nelson [1], is recently a subject of intensive theoretical and 
experimental research [2,3]. Such a kaon-condensed phase of dense matter is 
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of direct astrophysical relevance as it would form in neutron stars and could 
affect their properties. In this paper we study how sensitive is the cooling of 
neutron stars to the presence of the kaon condensates. A particularly impor- 
tant problem in this regard, which we address here, is to assess the influence 
of the uncertainty of the high density behaviour of nuclear symmetry energy 
on the charged kaon condensation itself. Other strange hadrons, like neutral 
kaons and hyperons could also appear in neutron star matter. Their inclusion 
here, however, would blur the role of the nuclear symmetry energy as it would 
introduce additional uncertainty due to poorly known nucleon-hyperon inter- 
actions. The analysis reported here will be extended to account for hyperons 
in a forthcoming paper. Model calculations show that hyperons are likely to 
influence the results concerning kaon condensation [4]. 

Among theoretical approaches we may distinguish those based on the chiral 
theory [1] and the ones using the meson exchange picture. The former ones 
have their roots in the spontaneously broken SU{3)xSU{3) symmetry. It gives 
a solid base to study kaon-nucleon interactions where kaons are treated as 
pseudo-Goldstone bosons. In the meson exchange picture, the kaon couplings 
to nucleons are realized indirectly through the exchange of vector and scalar 
mesons [5]. The predictive power of chiral theory is lost in such approaches. 
In order to study the kaon condensation in the neutron star matter a realistic 
model of nucleon-nucleon interactions should be considered. 

The nucleon-nucleon part of the interaction Hamiltonian is usually taken from 
other theories of nuclear matter, such as e.g. many-body models with realistic 
potentials or relativistic mean field (RMF) theory. Realistic potential approach 
is firmly based on experiment. Two-nucleon potentials are fitted to the scat- 
tering data. Three-body potentials must be added to get agreement with the 
binding energy of light nuclei ("^H, "^He) and the saturation properties of nu- 
clear matter. The energy of nuclear matter is usually derived in a variational 
approach [6] or by using other techniques as e.g. some version of the Brueckner 
theory. 

The RMF approach starts from postulating the Lagrangian for nucleons and 
mesons which takes into account relevant internal symmetries (as e.g. isospin). 
The meson-nucleon interactions are assumed to be of a Yukawa type with 
coupling constants treated as free parameters which are fitted to the saturation 
point properties obtained from the mass formula and from breathing modes 
of giant resonances [7] . 

Unfortunately, up to now there is no consistent description of the ground 
state properties of nuclear matter in terms of the chiral theory with its free- 
space parameters. Recently some novel approaches were proposed. One of them 
starts from the chiral Lagrangian [8] but still treats its crucial parameters in 
the spirit of RMF model - fitting them to the saturation point properties. A 
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different idea l:)ased on the power counting scheme was presented in [9], but 
its extrapolation to higher densities seems not to be well justified. 

In the next section we discuss the high density behaviour of the nuclear symme- 
try energy. One should stress that the symmetry energy governs the chemical 

composition of nuclear matter in the beta equilibrium and hence it determines 
the particle species present in neutron stars. In Sect. 3 the chiral model kaon- 
nucleon interactions is derived. In Sect. 4 the kaon condensate properties and 
the role of nuclear symmetry energy in a kaon-condensed matter are studied. 
In particular we assess how uncertainty of the behaviour of the symmetry 
energy at high densities influences conclusions regarding the formation and 
properties of the kaon-condensed neutron star matter. Implications for the 
cooling of neutron stars are investigated in Sect. 5. 



2 Nuclear symmetry energy 

A common feature of the RMF models is that the nuclear symmetry energy, 
whose value at the saturation density is well known to be about 30 MeV, 
strongly increases with density of the nuclear matter. It is not the case for 
realistic potential models of nuclear matter in the variational approach [6], 
where the symmetry energy saturates and then bends over at higher densities. 
The high-density behaviour of symmetry energy is the worst known property 
of dense matter [10,12], with different nuclear models giving contradictory 
predictions, as in the case of the RMF theory and realistic potential models 
which provide incompatible predictions of the nuclear symmetry energy at high 
densities. Such divergent predictions are the source of a major uncertainty of 
high density properties of neutron star matter [11], including the formation of 
kaon condensates [16]. 

To assess the influence of the symmetry energy we explore various forms of its 
high density behaviour as provided by different models of dense nuclear matter. 
We use three types of the high density symmetry energy obtained from the 
realistic potential theory [6] and also the symmetry energy derived from the 
RMF model which are displayed in Fig.l. Some preliminary results are given in 
[16], where it was shown that low values of the symmetry energy may exclude 
the condensation. However, if the condensate is formed, the proton fraction 
goes very high, even to complete protonization of nuclear matter, a; = 1. These 
results encourage us to explore the implications for cooling of neutron stars. 

The nuclear symmetry energy is defined as 
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where E{n, x) is the energy per baryon as a function of the baryon number 
density, n, and the proton fraction, x. The nuclear symmetry energy plays an 
essential role in the neutron star matter when the beta equilibrium among 
all its components is achieved. Prom thermodynamic identities it follows that 
dE/dx — /In — l-ip, where /in, IJ'p are neutron and proton chemical potentials. 
The isospin symmetry allows us to expand E{n, x) in even powers of (x — |). 
The expansion up to the second order is enough, i.e. the term ~ (x — seems 
to be negligible [13]. We may thus write with a sufficient accuracy 

4(1 - 2x)E, = - i^p. (1) 



The beta equilibrium yields a sequence of equalities between chemical po- 
tentials for particles present in the system and hence the symmetry energy 
governs the chemical composition of the neutron star matter. For example, in 
the simplest case of n,p, e matter, we have 4(1 — 2x)Es = fin — fJ^p = A*e- Tak- 
ing into account the neutrality condition, kg = kp, and neglecting the electron 
mass we obtain the proton fraction 

^^1(1__^). (2) 



At the saturation density, where the value of symmetry energy is known to be 
Es{no) = 04 = 30 MeV [14], the proton abundance is about 4%. As the proton 
Fermi momentum is slowly varying with density, kp ~ n^^^, we conclude that if 
Eg increases faster than n^/^ then the neutron star matter becomes more and 
more symmetric at high densities: a; |, (muons, produced when k^ > m^, 
make this effect even stronger because protons have to neutralize additional 
negative charge carriers). This is the case of RMF models. In its basic version 
where the only isovector interaction is due to Yukawa coupling to the p meson, 
the symmetry energy takes the form [15]: 

EJn) = + jLn, kF^i^nfl^ (3) 



The first term is purely kinematic, coming from the energy difference of proton 
and neutron Fermi seas. The second term comes just from interaction and is 
determined by the coupling constant. As we see. Eg grows slightly faster than 
linear with density and this is its typical behaviour in RMF theories ^ . In an 
opposite case, when the symmetry energy decreases with density, one expects 
the neutron star matter to be more neutron-rich, a; — > 0, at high densities. 
This simple picture need not be always true if kaons appear. Although kaons, 

^ In our earlier work [15] we tested if the inclusion of isovector-scalar component 
due to the 5 meson, could change the rapid growth of the symmetry energy. 



4 



similarly as electrons and muons, carry negative charge, their presence in the 
neutron star matter may have unexpected consequences. 



3 Model for K — N interactions 



In our investigation we use the effective chiral theory proposed by Kaplan and 
Nelson [1] and extensively explored in subsequent works. The Lagrangian for 
kaons represents a nonlinear effective theory which is not renormalizable, but 
it makes a good starting point for chiral perturbation theory [17] (we use the 
same symbols as in [18], for details see Appendix A.l). 



^Trd^Ud^'U^ + TrB{iYDf, - mB)B 
+ F TrBYldA, B]+ D TrB^^'^.iA,, B} 
+c{VrM{U + U^) + a{TrB{^M^ + ^'^MS})B 
+a2TrBB{iMi + ^'^ M^) + a^^IrBB 1r{^ME, 



(4) 



The Lagrangian consists of one part which is SU{3) x SU{3) symmetric (the 
first four terms in (4)), and the other part which breaks the chiral symmetry, 
where the breaking strength is determined by the parameters ai, a2, 03, and c, 
as well as by the quark mass matrix M. which is 



'^O ^ 



M 



00 

m, 



(5) 



i.e. the light quarks are assumed to be massless. Parameters in the symmetric 
part are / = 93 MeV, D = 0.81, F = 0.44. In the symmetry-breaking sector 
c represents the vacuum expectation value of the scalar quark density qq [19]. 
The Gell-Mann-Oakes-Renner relation states that mj^ — 2cms/p. Mass split- 
tings for the strange baryons determine the values of aim^, 02/^5 and we use 
after [20] airris = —67 MeV and a2ms = 134 MeV. The last parameter a^rng 
is subject to a large uncertainty. In principle it should be extracted from the 
kaon-nucleon scattering or kaonic atoms data [2] through the kaon-nucleon 
sigma term T^kn — ^ms{N\uu + ss\N) = — (|ai + a2 + 2a^ms , but up to 
now its value has not been determined. To assess its influence on further re- 
sults we adopt its values in the range —310 MeV < flsm^ < — 134 MeV which 
corresponds to a reasonable strangeness content of the proton in the range: 
0.2 > {ss)p/ {uu + dd + ss)p > [21]. The value of a^rUg appears to be impor- 
tant for the kaon condensation because the kaon-nucleon sigma term provides 
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the attractive interaction between kaons and nucleons, and makes the kaon 
effective mass smaller in the nuclear medium. Negatively charged kaons con- 
dense in the nucleon matter if their effective mass becomes smaller then the 
electron chemical potential /ig. Assuming that the only relevant kaon- nucleon 
interaction is in the s-wave, the expectation value of the kaon field {K~) is 
spatially uniform. According to the Baym theorem [22], its time dependence 
is 

{R-} = ^exp{-^Kt), (6) 



where 6 is an unknown amplitude of the condensate. 

Formally, pure nucleon-nucleon interactions should also emerge from the chiral 
theory as it is supposed to represent the low energy limit of QCD. However, 
it was not derived up to now, so we use the chiral Lagrangian only to extract 
the kaon-nucleon part of the interactions and we adopt the nucleon-nucleon 
part from the realistic potential models of nuclear matter. This means that 
the energy density of our problem is a sum of three contributions: the kaon- 
nucleon interaction energy density (derived from restricted only to kaon 
and nucleon fields), the nucleon-nucleon energy density taken from [6] and the 
energy density for leptons (e, /x): 

The nucleon-nucleon contribution is 

Snn ^ -Ef'noi — ] +V{n)+mn + n{l-2xfEs{n), (8) 
5 Vno/ 



where 



V{n)=nV^{n), (9) 

Es{n)^\Ef {^^J'\VM- (10) 

The functions Vq and V2 depend on the parameterization of the nucleon- 
nucleon potential used to fit the experimental data. We made use of the func- 
tions given in [6] corresponding to UV14+UVII, AV14+UVII, UV14+TNI in- 
teractions, where UV14 and AV14 are 14-operators two-body potentials from 
Urbana and Argonne. The three-body potentials UVIII or TNI are important 
to reproduce the saturation point properties by providing a proper contribu- 
tion in the isoscalar sector. The UV14-I-TNI potentials reproduce best the 
saturation density no = 0.16 fm^ and the binding energy wq — MeV. One 
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Fig. 1. The nuclear symmetry energy for different realistic potential models: 
UV14+UVII, AV14+UVII, and UV14+TNL The line labelled Ei corresponds to 
the symmetry energy for the RMF model [18]. 

can notice that the three-body forces give different energies also in the isovec- 
tor sector, i.e. the symmetry energy, as shown in Fig.l. For the same two-body 
potential, UV14, the three-body potentials UVII and TNI give different ener- 
gies at higher densities, although both of them produce the same symmetry 
energy at saturation point Es{nQ) ^ 30 MeV. In this figure, also the symmetry 
energy, typical for RMF theory, is shown. It is given by the formula [18] 



where u — n/rio and the interaction contribution F{u) depends linearly on 
density, F{u) — u. The nuclear symmetry energies from reahstic potential 
models differs from one another but only quantitatively, whereas the RMF 
prediction Ei is qualitatively different - it never saturates at higher densities. 

The energy density for leptons is determined in a standard way through their 
Fermi momenta 



where the function g {t) comes from integration over the Fermi sea of massive 
fermions: 



o 



(11) 




(12) 



git) 



1 



{{2t^ + t)VlT^ - arsinhi) , 



(13) 
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Fig. 2. Energy per particle for kaon-condcnsed neutron star matter with UV14+TNI 
(upper panel) and AV14+UVII (lower panel) interactions. The labels a, b and c 
correspond to the values of asrug from Table 1. 




a; ; a; > 
; X < 



To derive the kaon-nucleon contribution, Ekn, in (7) we use the standard 
prescription to get the energy density from a given Lagrangian (Appendix 
A.l). Employing the Baym theorem we obtain finally 

9 9 
^KN = sin^ 9 + 2m]^p sin^ - + n{2aix + 2a2 + 403)771^ sin^ - . (14) 



4 Thermodynamics of the system 

At zero temperature the energy density should be canonically a function of 
the number densities e (Tip, 77.„, 77./^) (leptons may be disregarded without loss 
of generality). Instead of number densities we have n, in the expressions 
(14), (8) entering £, which describe the thermodynamic state completely. The 
fourth, superfluous variable, ^, is a free microscopic parameter with respect to 



8 



which the energy is minimized when all densities are fixed, 



(15) 



This condition leads to the equation for Q (Appendix A. 2) which now may be 
expressed in terms of the equilibrium state determined by n, hk-, 

cos9 = J' ^ {jn]<:P + \n{2aiX-\-2a2-\-^a^)ms - \iiKn{l-\-x)\ . (16) 

The function 6{n, x, fix) plays the role of the order parameter in the analysis 
of the phase transition. The kaon condensate develops provided the above 
equation possesses a solution. If there is no solution, i.e. the right-hand side 
of eq.(16) does not fall in the range (-1,1), this means dense matter is allowed 
to contain only the normal nucleon phase. In the neutron star matter the free 
thermodynamic parameters n,x, /ik, fJ'e, fJ'n a-re subject to further constraints 
which come from the beta equilibrium and charge neutrality. 

On the timescale relevant to cooling of neutron stars, the nucleon matter 
with a kaon condensate approaches the beta equilibrium through the following 
reactions : 

n-^p + l + vi 

n^p + K~ , (17) 
/ K~ + 1/1 , I — e, II 

As neutrinos are supposed to leave freely, from the above reactions one derives 
corresponding relations among chemical potentials for all particles present in 
the system: 

I^K ^ IJ-e ^ IJ-ix ^ IJ-n - IJ-p (18) 

The above sequence of equalities shows that there are two independent chem- 
ical potentials, let's say and fi^- These chemical potentials represent two 
quantities, the electric charge and the baryon number, which are conserved in 
the beta- reactions (17). This means that the charge and the baryon number 
uniquely describe the thermodynamic state of our system. In further calcula- 
tions we put 

11 = He , 

and 11 will represent the chemical potential of all negatively charged particles. 
Chemical potentials for leptons are determined by their Fermi momenta, /ig ~ 
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ke, and, fj,^ — ^k'^ + m^, and for nucleons they must be derived through a 

proper differentiation, fjip^n — {de/dnp^n)ni, i-e. with all remaining densities 
fixed: 



I_i^eP (—y^\m + V'{n) + n{l-2xfE'^{n) + {l-Ax^)E,{n) 
Vno/ 



+ {2a2ms + AasTHs - /ik) sin^ ^, (19) 

Hj^E'^f' — + m + V'{n) + n{l-2x)^E'^{n) - {l-2x){3-2x)Es{n) 

+ {2aims + 2a2ms + 403771.5 — 2iik) sin^ -■ (20) 

With these expressions for //p, //p, the beta equilibrium condition (18) leads to 
the equation 



/X = 4(1 - 2x)Es{n) I cos^ - - 2aim, tan^ -. (21) 

Zj Zi 



While the beta equihbrium condition is easily implemented, the assumption as 
to the charge neutrality needs further qualification. We must make distinction 
between local and global electric neutrahty of the system. For locally neutral 
matter in any point densities of charged particles satisfy the condition 

rip-riK -ne-ni^ = 0. 

The only unknown quantity in the above formula is the negative kaon density 
TiK, which is obtained from the conserved electric current corresponding to 
the U{1) symmetry of the kaon-nucleon Lagrangian 

riK = i{pK+K^ -pK-K~), 



where Pk+ , Pk- are the canonical momenta. After somewhat laborious calcu- 
lations they may be derived directly from the Lagrangian (Appendix A.l). 
Making use of the Baym theorem we get the final expression for the charge 
density carried by kaons, 



riK^ f pi 8111^9 + n{l + x)sm'^- , (22) 

z 



and the local charge neutrality condition becomes: 



= /Vsin^ e + n(l + x) sm' ^ + ^ + 0(l^)^-^^-^^ - nx . (23) 
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Fig. 3. The kaon chemical potential for indicated nuclear models and different values 
of asms (models a,b and c from Table 1). 

When a phase transition takes place, the charge separation in two different 
phases preserving merely the neutrality of the matter as a whole, but not 
on a small scale, can occur, as it was shown in [23]. Then the equation (23) 
does not longer hold. The problem of mixed phase formation will be discussed 
elsewhere. 



The equations (23,21,16) allow us to find the total energy density (7) as a 
function of the baryon density only, e{n). The onset of condensation occurs 
at a critical density ric where the function 6 starts to deviate from zero and 
the energy density of the phase with the condensate becomes smaller than 
that of normal matter, Fig.2. As it was shown in [16], the critical density, 
ric, depends sensitively on the shape of the symmetry energy, Es{n), with the 
eventual exclusion of condensation if the strangeness is too small. The values 
of ric for different models of the symmetry energy. Eg, and some representative 
values of strangeness are presented in Table 1. Likewise, the other quantities. 

Table 1 

Critical density for the kaon condensation in units of tiq. In case of UV14+TNI 
interactions with the minimal strangeness 037^5 = — 134MeV the kaon condensate 
does not appear. 





asms [MeV] 


El 


UV14+UVII 


AV14-FUVII 


UV14-FTNI 


a 


-310 


2.4 


2.6 


2.7 


2.8 


b 


-230 


3.0 


3.4 


3.7 


4.1 


c 


-134 


4.2 


5.0 


5.9 





important for the structure of neutron stars, change radically with Es- Para- 
doxically, low values of the symmetry energy (which usually reduce x) in the 
presence of the condensation lift the proton abundance even up to a complete 
protonization, x — 1. In the most extremal case of the UV14-I-TNI interac- 
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tions the protonization sets in just after formation of the condensate. This is 
because kaons appear at a small value of the kaon chemical potential fi, that 
means they may be easily produced as the chemical potential measures the 
cost of energy to add one particle to the system. The kaon chemical potential 
in the condensate is shown in Fig. 3 for some models we study here. 



With the presence of kaon condensate the overall chemical composition of the 
neutron star matter becomes somewhat exotic. This is shown in Fig.4 and 5. 
Let us consider first the case of AV14+UVII (UV14+UVII gives similar results 
but less pronounced). Kaons, as bosons, are produced in the zero momentum 
state unlike leptons which must be put on the top of the Fermi sea. The 
negative charge carried by kaons is so abundant that the kaons supersede 
electrons and muons in dense matter - the curves for x and xk overlap (case (c) 
in Fig.4). When strangeness is high (case (a) in Fig.4) the kaons are produced 
at lower density, xk approaches the value xk — 1 and even may exceed unity. 
If Xk is so high, positive leptons {xe+^ < 0) have to appear to ensure charge 
neutrality. With further increase of density the situation changes again. As 
the proton fraction x cannot be greater than 1, the equation (16) is no longer 
valid. Strictly speaking, the energy has no minimum for x inside the allowed 
range, x e (0, 1). So, it takes its minimal value at the edge of the range of x 
which is just at x — 1. Then stops to decrease, positive leptons disappear 
and negative leptons may be restored. This fact is better seen in Fig. 5. In 
the UV14+TNI case positive leptons appear as well but with much smaller 
amount because the effect of a; = 1 constraint is very strong. And then, at 
higher density, the composition of matter does not depend on the strangeness 
value. The curves for xk and Xg+^u overlap for high n (see Fig. 5). 



The above discussion clearly indicates the situation becomes complex. There 
is a common conviction that high proton fraction a; > 0.15 ensures very fast 
coohng of neutron stars by direct URCA processes [24]. This statement is 
true only for usual n, p, e, /i matter. Direct URCA cycles take place if Fermi 
momenta for protons, neutrons and leptons form a closed triangle. All the 
Fermi momenta may be expressed in terms of x. So, the triangle condition leads 
to the existence of the threshold value, a;uRCA, which lies between 0.11 and 
0.15, depending on the presence of muons. For nuclear matter with kaons we 
have an additional degree of freedom - 9, which determines the kaon abundance 
altogether with fi (22). This is the reason why we cannot formulate any simple 
answer to the question if direct URCA processes are allowed above some value 
of X or not. In our case, the triangle condition must be examined explicitly in 
the whole range of densities. 
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Fig. 4. The composition of the kaon-condensed neutron star matter for AV14+UVII 
potentials. Two families of curves indicated as (a) and (c) correspond to maximal 
and zero strangeness in Table 1, respectively. Bold lines represent the proton fraction 
X, thin lines xk and .Xe+^. For comparison, dashed lines correspond the proton 
fraction x calculated with the same function Vq and the Ei model of the nuclear 
symmetry symmetry energy. 




Fig. 5. The composition of the kaon-condensed neutron star matter for UV14+TNI 
potentials. Line labels as in Fig. 4. For zero strangeness (c) there is no condensation, 

so only (a) and (b) families are presented. At high density all particle fractions Xi 
are independent of the strangeness and corresponding curves overlap for different 
families. 

5 Cooling of matter with kaon condensate 



When a kaon condensate is formed, besides the usual direct beta reactions 
there are also possible reactions which are specific to the kaon-condensed neu- 
tron star matter. The direct URCA cycles take place among quasi-particles 
p,h 
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i (24) 
p + I ^ h + vi, 

and we call them dURCA. The states \n) are linear combinations of pure 
proton and neutron state which depend on the condensate amplitude 6 [25]. 
As the quasi-particles have no definite charge one may write another admitted 
type of reactions 



n^h + l+ ui, (25) 
p^p + l+ vi . (26) 

Microscopically, it means that the necessary charge is taken from the conden- 
sate, e.g. n + K ^ n + l + vi and thus they are called kaon induced UKC A. For 
all particles taking part in the above reactions, the momentum conservation 
gives the following relations among their Fermi momenta: 



kURCA : 2kn > h or 2kp > k (27) 
dURCA : \kp- kn\ < ki < kp + kn ■ (28) 

The neutrino momentum, which is k,y ~ kT < IMeV, may be neglected as it 
is much smaller than momenta of the other particles. Likewise, kaons do not 
enter into the relation (28) as the condensation takes place in the s-wave and 
the kaon momenta do not contribute. Let's focus on the direct URCA cycle. 
The three terms of the double inequality are presented for various models in 
Fig. 6. When negative leptons are present, the right-hand side inequality in 
(28) is always fulfilled , because the local neutrality means that k^ < kp. It 
does not have to be true when positive leptons appear. In principle, although 
very unlikely, it is possible that kg > kp and for clarity we put the sum kp ~\~ kfi 
on the plots as well. The main conclusion from the calculations shown in Fig. 6 
is that the exceeding by the proton fraction the threshold value xurca does 
not ensure the presence of the dURCA cycle. Too high values of x block this 
cycle due to small fraction of neutrons. Also, a non-typical behaviour of the 
chemical potential /i, which crosses zero at some density (Fig. 3), suppresses 
dURCA in a narrow zone due to small lepton abundance. From the plots in 
Fig. 6 the tendency of suppressing the dURCA process with smaller values of 
the symmetry energy can be inferred. However, the importance of this effect 
is not so great, because always one of the kaon-induced processes (with a pro- 
ton or a neutron) operates as may be easily concluded from Fig.6. The sum 
kp + kfi is always greater than kg so, either 2kp > kg or 2kn > kg or both of 
them are fulfilled. In the thermal evolution of neutron stars we distinguish 
two cooling scenarios: fast - driven by direct URCA and slow - driven by 
modified URCA processes. Modified URCA always operate, as they are not 
limited by any kinematical constraints, however their emissivities are much 
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Fig. 6. The dURCA triangle condition for various nuclear models and different values 

of strangeness. The solid lines correspond to kp + kn, dashed lines to kg and dotted 
lines to \kp — kn\. The shaded areas represent regions where dURCA is allowed 
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Table 2 

Emissivities of the three fermion processes in the kaon-condensed neutron star mat- 
ter. 



cycle 


reaction 


hi -fuRCA 


dURCA 


h—>-p + l + ui 


cos^ (^/2) 


kURCA 


h ^ h + l + ui 
■p^p + l + vi 


\ sin^e" tan^ 
sin^6' tan^ 6'c 




Fig. 7. The total neutrino emissivity in URCA processes for various nuclear models. 



smaller, about 5 to 6 orders of magnitude less than the direct ones. So they 
provide much slower cooling. The kaon-induced URCA processes (like direct 
URCA) belong to the class of three-fermion processes but their emissivity 
is smaller because they proceed through strangeness-changing reactions. The 
rate of such reactions scales as sin^c, whereas direct URCA scales as cos^c, 
where is the Cabibbo angle. Both kinds of processes depend on the kaon 
abundance which is described by the condensate amplitude Q [26]. In Table 2 
the ratio of emissivities of kaon-condensed neutron star matter to emissivity 
by direct URCA of normal matter with no condensate is presented. The emis- 
sivity finally depends on the baryon number density through /urca and the 
condensate amplitude Q. The total emissivity is a sum of all the contributions 
from direct and kaon-induced URCA cycles with electrons as well as muons 
which are present if //e > In Fig. 7 the total emissivity is shown. 
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6 Summciry and discussion 



We have shown that the formation and properties of the kaon-condensed phase 
of neutron star matter are quite sensitive to the high density behaviour of the 
nuclear symmetry energy Es{n), which still remains the most poorly known 
property of dense matter. In particular, for Es decreasing to negative values the 
formation of the kaon condensate can be inhibited for the lowest absolute value 
of asTTis. This is in stark contrast to the results for monotonically increasing 
Es{n), when the condensation occurs for any value of Osm^. The properties of 
the kaon-condensed neutron star matter and the abundance of particle species 
are sensitive to the high density form of the nuclear symmetry energy as shown 
in Figs. 4-5. A surprising finding is that the direct URCA process is often not 
allowed in the kaon-condensed matter for decreasing Es{n). 

The above results allow us to assess how great is the influence of nuclear 
force models on the total neutrino emissivity of kaon-condensed neutron star 
matter. In case of increasing symmetry energy as in the RMF model (Ei) 
dURCA operates almost always when kaon condensate is formed, except of 
a narrow zone around density where the chemical potential is zero. For 
realistic potential models dURCA is systematically suppressed to more narrow 
density range with decreasing symmetry energy. For an extreme case of the 
UV14-I-TNI interactions, the allowed zone shrinks eventually to a very narrow 
vicinity of some density close to the threshold value for condensation ric- The 
kaon-induced URCA branch is always present but at the level of about 10 
times smaller, as shown in Fig. 7. 

Another interesting question concerns the role of As was shown in [27] neu- 
tral kaons arc easily produced just after the onset of K~ condensation. Then 
the kaon-condensed matter becomes even more effectively isospin-symmetrized 
than in the absence of K^. We would like to stress that this effect is sensitive 
to the form of the nuclear symmetry energy. The effect occurring for increas- 
ing Es{n) is not observed in the case of decreasing Es{n). High values of x 
which we obtain in our model make the condensation more difficult. The 
dispersion relations for anti-kaons are: 



n{l + x) If n{l + x)Y , 2 , /o , n , A \ ^ 



1/2 



= V/T^ + y[ 4J2 j +^K+ (2aix + 2a2 + 4a3)m,— | (29 

n(2-x) ((n{2-x)Y o ^ ^ 

= \.o + \ro +rn\ + (2ai(l - x) + 2a2 + 4a3)m, 



4y2 \\ 4/^ / -K'v j^v /' ^ • o/*2j2 



The first term in the expression (30) for uj^o, which is the leading term, in- 
cludes positive contribution from x, so high values of x make the slope of Wj^o, 
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as a function of n, more fiat after tlie cliarged iiaons appear. Fig.8 presents tliis 
effect ciearly. We ieave more detaiied discussion of the neutral kaon condensate 
for future work. 

AV14+UVII a3ms=-134 MeV 

1 

0.8 
^0.6 
0.4 
0.2 



2 4 6 8 10 




2 4 6 8 10 

n/no 



Fig. 8. The proton fraction x and anti-kaons energy in dense matter. The leveling 
off of (jJk for corresponds to the onset of K~ condensation which drives the 
strong increase of the proton fraction shown in the upper panel. 

Note added: After completion of this work dr. T. Muto brought to our at- 
tention his research [28] concernig kaon condensation and cooling of neutron 
stars. 




18 



A Calculation of thermodynamic quantities 



A . 1 Basic formulae 



/-M 

A matrix U e S'[/(3) is defined asU — exp(iv2— ), where M represents the 
pseudoscalar meson octet 









7r+ 




















^0 





(A.l) 



Matrix ^ appearing in (4) is such that — U and B corresponds to the baryon 
octet 



B - 



V 



■^0 



-VIA/ 



(A.2) 



When the only non-zero elements of the meson matrix M are K'^ and K 
fields, then U has the form (below we use dimensionless meson fields A;^ = 



I 



u 



COS V k k+ 




ifc+ sin 



Vk-k+ 

1 



(A.3) 



Hence the canonical momenta corresponding to the rescaled kaon fields p± 
dC 



sm 



2 X 



P+ = -i (nn + 2pp)——^k 



sm 



k+ 

2 X 



(A.4) 



P-— i{nn + 2pp) :^k~^ 

2x^ 
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f ( k+dtk- + k-dtk+ k-dtk+ - k+dtk' sin^x 
A [ k- k- . 



where x = Vk+k~. The energy density for the kaon-nucleon sector is given by 
the Hamiltonian Hkn — P+dtk'^ — p-dtk" — and it reads 



^ f(k+d,k- + k-d,k* f p(k+d,k- -k-d,k+ f sin' X ,,„, 

n,.= ^ (A.6) 

+ Ozf^mj^ + (2a2 + 403) rris nn + (2ai + 2a2 + 403) nispp] sin^ — 
The charge density uk — i{p-\~k'^ —p^k~) becomes 



UK = {nn + 2pp) sin2 | + '4-{k+d^k- - k'd^k+f-^ (A.7) 

^ ^ X 



The Baym theorem for k"^ fields allows us to write them as = ^ exp='=*'*^* 
and with this form we get the equation (14,22). 

A. 2 Minimum oj energy 

The derivative of Si^Arwith respect to Q, keeping = n^x^nx fixed, is 
(deKN\ deKN{n,x,iiK,0) , deKN{n,x,iiK,0) fdfj,K\ , . 

[^)„r + a^, i^j„; 

It remains to find ^ qq ^ ' '^^^^ ^^^^ ^® done by taking derivative of both 
sides of (22) 

= /^sin2^(^) + 2/Vi^sin^cos^ + in(l + x) sin ^ . 



Hi 



This approach is, of course, equivalent to minimization of a new thermody- 
namic potential e = Skn — /J-nx with respect to 9 with n,x,ii fixed, as was 
done in [18]. 
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